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ABSTRACT 

The energy requirements of gamma ray bursts have in past been poorly constrained because of three 
major uncertainties: The distances to bursts, the degree of burst beaming, and the efficiency of gamma 
ray production. The first of these has been resolved, with both indirect evidence (the distribution of 
bursts in flux and position) and direct evidence (redshifted absorption features in the afterglow spectrum 
of GRB 970508) pointing to cosmological distances. We now wish to address the second uncertainty. 
Afterglows allow a statistical test of beaming, described in an earlier paper. In this paper, we modify 
a standard fireball afterglow model to explore the effects of beaming on burst remnant dynamics and 
afterglow emission. If the burst ejecta are beamed into angle £ m , the burst remnant's evolution changes 
qualitatively once its bulk Lorentz factor T <J l/£m : Before this, T declines as a power law of radius, 
while afterwards, it declines exponentially. This change results in a broken power law light curve whose 
late-time decay is faster than expected for a purely spherical geometry. These predictions disagree with 
afterglow observations of GRB 970508. We explored several variations on our model, but none seems 
able to change this result. We therefore suggest that this burst is unlikely to have been highly beamed, 
and that its energy requirements were near those of isotropic models. More recent afterglows may offer 
the first practical applications for our beamed models. 



1. INTRODUCTION 

Understanding the energy requirements and event rates 
of gamma ray bursts is necessary for any quantitative eval- 
uation of a candidate burst progenitor. We need to know 
both how many progenitors we expect, and how much en- 
ergy they need to produce in a single event. Until recently, 
both quantities were uncertain to ~ 10 orders of magni- 
tude because of the unknown distance to the bursts. The 
afterglow of GRB 970508 effectively ended that debate, be- 
cause it showed absorption lines at a cosmological redshift 
= 0.835; Metzger et al 1997). This builds on earlier 
results from the Burst and Transient Source Experiment 
(BATSE), which showed that the burst distribution on the 
sky is exquisitely isotropic while the distribution in flux is 
inhomogeneous (Meegan et al 1996). These observations 
are best explained if the bursts are at cosmological dis- 
tances. A very extended Galactic halo distribution might 
also work, but it would have to be unlike any other known 
population of Galactic objects. The isotropy is perhaps 
most important now for showing that multiple-population 
scenarios for gamma ray bursts cannot put any substantial 
fraction of the bursters at Galactic distances. It thus con- 
nects the GRB 970508 redshift bound to the vast majority 
of the burst population. 

The dominant remaining uncertainty in the bursters' 
energy requirements is now whether the bursts radiate 
isotropically or are beamed into a very small solid an- 
gle. Such beaming is allowed (though not required) by the 
gamma ray observations, because the ejecta from gamma 
ray bursts must be highly relativistic to explain the spec- 
tral properties of the emergent radiation (Paczyhski 1986, 



Goodman 1986), with inferred minimum Lorentz factors 
T ^ 100 (Woods & Loeb 1995). The gamma rays we ob- 
serve are therefore only those from material moving within 
angle 1/r of the line of sight, and offer no straightforward 
way of determining whether there are eject outside this 
narrow cone. 

These large Lorentz factors lead naturally to predictions 
of afterglow emission at longer wavelengths as the burst 
ejecta decelerate and interact with the surrounding mate- 
rial (Paczyhski & Rhoads 1993; Katz 1994; Meszaros & 
Rees 1997a). The characteristic frequency for this after- 
glow emission depends on the Lorentz factor of the burst 
remnant, and both decrease as the remnant evolves. Such 
models scored a recent triumph with the detection of X- 
ray, optical, and radio afterglows from gamma ray bursts 
(GRBs) early in 1997 (e.g., Costa et al 1997; van Paradijs 
et al 1997; Bond 1997; Frail et al 1997). The observed 
properties of the transients are in good overall agreement 
with the predictions of afterglow models (Wijers, Rees, & 
Meszaros 1997; Waxman 1997a, b), although some worries 
remain (Dar 1997). 

Because beaming depends on the relativistic nature of 
the flow, afterglows can be used to test the burst beaming 
hypothesis. At least two such tests are possible. First, 
because T is lower at the time of afterglow emission than 
during the GRB itself, the afterglow cannot be as colli- 
mated as the GRB can. This implies that the afterglow 
event rate should exceed the GRB event rate substantially 
if bursts are strongly beamed. Allowing for finite detection 
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thresholds, 



where N\ , N% are the measured event rates above our de- 
tection thresholds at our two frequencies; N\2 is the rate of 
events above threshold at both frequencies; and Cli, f2 2 are 
the solid angles into which emission is beamed at the two 
frequencies. A full derivation of this result and discussion 
of its application is given in Rhoads (1997a). 

The second test is based on differences between the dy- 
namical evolution of beamed and isotropic bursts. Burst 
ejecta decelerate through their interaction with the am- 
bient medium. If the ejecta are initially beamed into a 
cone of opening angle ( m , the deceleration changes qual- 
itatively when the bulk Lorentz factor T drops to 1/Cm- 
Prior to this, the working surface (i.e. the area over which 
the expanding blast wave interacts with the surrounding 
medium) scales as r 2 . At later times, the ejecta cloud 
has undergone significant lateral expansion in its frame, 
and the working surface increases more rapidly with r, 
eventually approaching an exponential growth. Spherical 
symmetry prevents this transition from occurring in un- 
beamed bursts. A brief analysis of this effect was presented 
in Rhoads (1997b). 

We have two major aims in this paper. First, we will 
present a full derivation of the late time burst remnant 
dynamics for a beamed gamma ray burst. We support 
this by calculating the emergent synchrotron radiation for 
two electron energy distribution models, but we do not at- 
tempt to do so for all possible fireball emission scenarios. 
Second, we observe that our model is not consistent with 
any small-angle beaming of GRB 970508. This implies a 
substantial minimum energy for this burst. If radiative 
efficiencies are lower than ~ 10%, this limit approaches 
the maximum energy available in compact object merger 
events. We explore possible ways to evade this minimum 
energy requirement through other forms of beaming mod- 
els, but find none. We therefore conjecture that such mod- 
els cannot be constructed for GRB 970508 unless the usual 
fireball model assumptions about rclativistic blast wave 
physics are substantially modified, and challenge the com- 
munity to prove this assertion right or wrong. 

We explore the dynamical evolution of a model beamed 
burst in section In section |^ we incorporate a model for 
the electron energy spectrum and magnetic field strength 
and so predict the emergent synchrotron radiation. In sec- 
tion ^, we compare the model with observed afterglows. 
The early (1997) data appeared inconsistent with the 
beaming model, suggesting that bursts are fairly isotr opic 
and therefore very energetic events. Finally, in section |4.l| , 
we explore variations on our model to try to reduce the in- 
ferred energy needs of GRB 970508. We comment briefly 
on more recent data and summarize our conclusions in 
section g. 

2. DYNAMICAL CONSEQUENCES OF BEAMING 

We explore the effects of beaming on burst evolution 
using the notation of Paczyhski & Rhoads (1993). Let Tq 
and Mo be the initial Lorentz factor and ejecta mass, and 
£ m the opening angle into which the ejecta move. The 
burst energy is Eq = TqMqc 2 . Let r be the radial coor- 
dinate in the burster frame; t, t co , and the time from 



the event measured in the burster frame, comoving ejecta 
frame, and terrestrial observer's frame; and / the ratio of 
swept up mass to Mq. 

The key assumptions in our beamed burst model are 
that (1) the initial energy and mass per unit solid angle 
are constant at angles < Cm from the jet axis and zero 
for 6 > ( m ; (2) the total energy in the ejecta + swept- 
up material is approximately conserved; (3) the ambient 
medium has uniform density p; and (4) the cloud of ejecta 
+ swept-up material expands in its comoving frame at the 
sound speed c s = c/\/3 appropriate for relativistic mat- 
ter. The last of these assumptions implies that the work- 
ing surface of the expanding remnant has a transverse size 
~ Cm r + c s t co . The evolution of the burst changes when 
the second term dominates over the first. 

Each of these assumptions may be varied, but we be- 
lieve the qualitative change in burst remnant evolution 
will remain over a wide range of possible beaming mod- 
els. Removing assumption (4) is the only obvious way to 
turn off the dynamical effects of beaming, and even then 
observable breaks in the light curve are expected when 

r ~ l/Cm- 

There are several models in the literature that use ra- 
diative rather than adiabatic models, dropping our second 
assumption. The case for radiative bursts depends on the 
efficiency with which relativistic shocks transfer bulk ki- 
netic energy to magnetic fields and electrons, and I regard 
the validity of assumption (2) as an open question. For a 
closer examination of this issue, I refer the reader to pa- 
pers by Vietri (1997a, b) and by Katz & Piran (1997a), who 
advocate radiative models; and to Waxman, Kulkarni, & 
Frail (1998), who defend the adiabatic model. Meszaros, 
Rees, & Wijers (1997) point out that the dynamical conse- 
quences (r oc r~ 3 ) of radiative models depend on equipar- 
tition between protons, electrons, and magnetic fields be- 
ing maintained at all times. Thus, a short electron cooling 
time will affect the afterglow radiation, but will not neces- 
sarily result in r cx r~ 3 . Sari (1997) considers corrections 
to the adiabatic burst evolution for modest energy losses. 

Models that do not use assumption (1) have been 
discussed by Meszaros, Rees, & Wijers (1997) and 
Panaitescu, Meszaros, & Rees (1998). Finally, assump- 
tion (3) has been dropped by several authors (Vietri 1997b; 
Meszaros, Rees, & Wijers 1997; Panaitescu, Meszaros, & 
Rees 1998) in favor of a more general power law density 
p cx r~ 9 . Such models complicate the beamed burst analy- 
sis and will change the form of T(r) relation but will leave 
intact the basic conclusion that T(r) changes qualitatively 
when T <; l/( m . 

2.1. Dynamical Calculations: Numerical Integrations 

Given these assumptions, the full equations describing 
the burst remnant's evolution are 

f =TT f r 2 n m (r)p(r)dr , (2) 

Mlo Jo 

^rn ~ 7T(Cm + C s t co /ct) 2 W 7t(Cto + t co /Vit) 2 , (3) 

r = (To + /) /Vi + 2r / + / 2 « v/ToAf , (4) 

t = r/c , t co = f dt'/T , 
Jo 
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and tr 



(1 



dt'/2T 2 



(5) 



Equation^ is derived in Paczyriski & Rhoads (1993) from 
conservation of energy and momentum, along with alge- 
braic simplifications of equations [| for the spherical case. 
The definition of i® here includes the cosmological time 
dilation factor (1 + z) for a source at redshift z. Equa- 
tion [3] is not strictly valid when £ m ^> 1, but we will accept 
this deficiency since the error thereby introduced is not a 
dominant uncertainty in our results. 

These equations can be solved by numerical integration 
to yield f(r), T(r), and i®(r). Figure 1 shows T(r) from 
such integrations for an illustrative pair of models (one 
beamed, one isotropic). 

2.2. Dynamical Calculations: Analytic Integrations 

The most interesting dynamical change introduced by 
beaming is a transition from a power law T oc r -3 / 2 to an 
exponentially decaying regime T oc exp(— r/r T ). We will 
first give a derivation of the power law behavior. 



2.2.1. Power Law Regime 

Consider the approximate evolution equations for the 
regime where (a) l/r < / < r , so that T ss y / T /2f; 
and (b) c s t co <; ( m r (corresponding to / <; 9r C„). 



7f 2 

df/dr w —{Cmr) p , dt co /drz 

M V c 2 T 



2/ 



dt®/dr w (1 + z) 



J_ 

cTn 



(6) 



The initial conditions are / = and t = t co = t® = at 
r = 0. So we can easily integrate and obtain 



, = <m P r 3 



3 Mo 



3M r 



1/2 



3E 



2ir(? n c 2 r 3 p 



1/2 



_ ( %< 2 rnP \ 1/2 r „ 2 _ 2 r 



V 75£ / 5 c r 

i (l + z )^£^ r 4 

® y ' 12 E 



whence 



r = 2~ 5 / 4 



3E \ l7 Vl + z x3/S 



Cta 



(7) 

(8) 

(9) 
(10) 

(11) 



By making the substitutions — > An and (1 + z) — > 1 
in these results, we recover the evolution of a spherically 
symmetric burst remnant derived by Paczyriski & Rhoads 
(1993). 



2.2.2. Exponential Regime 

To demonstrate the exponential behavior, consider the 
approximate evolution equations for the regime where (a) 
l/ r o < / < r , so that T w y/T /2f; and (b) c s t co > ( m r 
(corresponding to / ;> 9T Cm) : 



df/dr 
eft® jdr 



TT 

Mo 



2.2 

C s tcoP 



dt co /dr ; 



2/ 

c 2 r n 



(i + «) 



cTn 



(12) 



By forming the ratio (df / dr) / (dt co / dr) and isolating terms 
with / and with t co , it follows that 



y/fdf 



TT CC^VTo 



x t 2 co dt co 



c 3 P\/To" 



V2 Mo c ° 3V2 

This is easily integrated to obtain 
' TTyT^C c 2 p 



Mn 



X tcc/dtco 



f 



3/2 



\/8Mo 



(4 - ci) 



(13) 



(14) 



where ci is a constant of integration. Using the initial con- 
ditions for the exponential regime derived below (eqn. ^ 
- |l^), one can show that the constant of integration is 
Ci = — 25 E (m/ (An pel), which becomes negligible once 
c s t C o 3> Cmf- Equation FF3 then becomes / oc t 2 a , and we 
see from equations 12] and ^ that /, T, t co , and i® will all 
behave exponentially with r in this regime. Retaining the 
constants of proportionality, we find 



/ oc exp(2r/r r ) where 

2 



r M " 


1/3 


[ E ° 1 


p 




irc 2 p 



1/3 



(15) 



Further algebra yields T oc exp(— r/r T ), t co oc \fj 

exp(r/r r ), and t® oc / oc exp(2r/r r ), so that T oc i® . 
Thus, while the evolution of T(r) changes from a power 
law to an exponential at T ~ 1/Cmj the evolution of t©(r) 
changes similarly. The net result is that r(i®) has a power 
law form in both regimes, but with a break in the slope 

from r oc i® 3 ^ 8 when T > l/( m to T oc ^®^ 2 when 

r<l/Cm; _ 

The initial conditions for the exponential regime are 
approximately set by inserting the transition condition 
c s t co = CmCt into the evolution equations for the power 
law regime. Denoting the values at this break with the 
subscript b, we have c s t co ^ = C, m ctb = Cmfb, which we 
combine with equation to obtain 



_ / 75Eo_ 

n ~Wpc 2 



1/3 



and J), ^ ( CF, 



The corresponding values for T, t®, and t co are 

r b = ^i— 

• 3 XV3 5 8/3 



(16) 



(17) 



J> 



(1 + *) F 
1 TT 



64 c s 



and t 



co,b 



\8irpcl 



(19) 
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The evolution in the exponential regime is then approx- 
imated by 



/ = fb x exp{2(r-r 6 )/r r } 

t® = te,b x exp{2(r-r 6 )/r r } 

r = T fc x exp{-(r - n)/r T } 

t C o = tco,b x cxp {(r - r b )/r r } 



(20) 



A thought experiment that will help understand the on- 
set of the exponential decay of T with radius is to consider 
the shape of a GRB remnant in a pressureless, uniform 
ambient medium at late times (after all motions have be- 
come nonrelativistic) . In the spherical case, the blast wave 
will leave behind a spherical cavity. In the beamed geom- 
etry, the cavity will be conical near the burster, but will 
change shape at the radius where the lateral expansion of 
the remnant becomes important. At this point, the cone 
flares, and the mass swept up per unit distance begins to 
grow faster than r 2 . This corresponds to the onset of the 
exponential T(r) regime. The cone continues to become 
rapidly wider until it reaches the radius where the rem- 
nant becomes nonrelativistic. The final cavity resembles 
the bell of a trumpet. It is unclear whether such remnants 
would survive long enough to be observed in a realistic 
interstellar medium. 



3. EMERGENT RADIATION 

The Lorentz factor T is not directly observable, and we 
ultimately want to predict observables like the frequency of 
peak emission v&. m , the flux density F v ^_ m at ^®, TO , and 
the angular size 9 of the afterglow. To do so, we need to 
introduce a model for the emission mechanism. We will re- 
strict our attention to synchrotron emission, which is the 
leading candidate for GRB afterglow emission. We first 
consider the case of optically thin emission with a steep 
electron energy spectrum. This emission model is used 
in many recent afterglow models (e.g. Waxman 1997a,b). 
We then repeat the calculation for the emission model of 
Paczyhski & Rhoads (1993). 

3.1. General equations: Optically thin case 

Our dynamical model for burst remnant evolution gives 
the volume V and internal energy density Ui of the ejecta 
as a function of expansion radius r. Detailed predictions 
of synchrotron emission require the magnetic field strength 
and the electron energy spectrum. We assume that the en- 
ergy density in magnetic fields and in relativistic electrons 
are fixed fractions £ B and £ e of the total internal energy 
density. The magnetic field strength B follows immedi- 
ately: 



B 



(21) 

(N.b., we use the notation of Paczyhski & Rhoads 1993. 
Some other authors have instead defined £ B in terms of 
the magnetic field strength, such that B cx £ B in their 
models; care must therefore be taken in comparing scaling 
laws under these alternative notations.) 

The electron energy spectrum requires additional as- 
sumptions. We first follow Waxman's (1997a,b) assump- 
tions, to facilitate comparison of our results for beamed 
bursts with his for unbeamed bursts. In the frame of 
the expanding blast wave, the swept-up ambient medium 



appears as a relativistic wind having Lorentz factor T. 
We assume that the electrons from the ambient medium 
have their direction of motion randomized in the blast 
wave frame. Moreover, they may achieve some degree 
of equipartition with the protons. The typical random 
motion Lorentz factor "f e for the swept-up electrons in 
the blast wave frame is then in the range r <J 7 e <J 
0.5(m p /m e )T. In terms of the energy density fraction in 
electrons, (-f e ) w ^ e (m p /m e )T. We further assume that 
the electrons in the original ejecta mass are not heated 
appreciably, so that the number of relativistic electrons 
is N e = fMo/{fi e m p ) (where u. e is the mean molecular 
weight per electron) rather than (1 + f)Mo/ (/j, e m p ). We 
take the electron energy £ to be distributed as a power law 
N(£) = £- p for £ min < £ < £ max , where N(£)d£ is the 
number of electrons with energies between £ and £ + d£. 
Finally, we assume that p > 2, so that the total electron 
energy J £ max £N(£)d£ is dominated by electrons with £ w 

£ min , and 7 e!pcak ~ (j e ) ~ £, e (m p /m e )T ss £ min /(m e c 2 ). 

The optical depth to synchrotron self-absorption 
is assumed to be small at the characteristic syn- 
chrotron frequency corresponding to £ m - m — 7e. P oak?7ieC 2 . 
In the comoving frame, this frequency is v co m = 
0.29 x 3/(47r)(sina) 7e 2 pcak ei?/(™ e c) 0.29 x 

(3/16)7 2 pcak e J B/(m e c) (Pacholczyk 1970, Rybicki & 
Lightman 1979), where the calculation of the mean pitch 
angle (sin a) = 7r/4 assumes an isotropic distribution of 
electron velocities and a tangled magnetic field. Wijers 
& Galama (1998) have integrated over the power law 
distribution of electron energies to show that the peak 
comoving frame frequency for a power law energy distri- 
bution becomes {v co , m ) £ = 3a; p /(47r) x ll_ pca ^cB / {m e c), 
where function of the power law index p, and 

where 0.64 ;> x p J> 0.45 for 2 < p < 3. Below this peak 
frequency, the flux density rises as v 1 / 3 , while at higher 
frequencies it falls as v~ a where a = (p — l)/2 (e.g., 
Rybicki & Lightman 1979, Pacholczyk 1970). Three ad- 
ditional breaks may occur, corresponding to the highest 
electron energy attained in the shock, the electron energy 
above which cooling is important, and the frequency where 
synchrotron self- absorption becomes important. We will 
comment on the cooling break below, and will ignore the 
other two breaks for the present. 

We first estimate the observer-frame frequency ^©, m at 
which the spectrum peaks. This is 



(l+/3cos(0))*r x (v co , m ) £ /(l + z) 



1 



1 + z 

1 Xp 
1 + Z 7T 



4F x 3*p x 7 e , P ea k e£ (22) 
3 47r m e c 



m p 



m e c 



where the factor (1 + /3 cos(#))r is the Lorentz transforma- 
tion for frequency, and where = \/l — F~ 2 is the expan- 
sion velocity as a fraction of lightspeed c. 9 is the angle 
between the velocity vector of radiating material and the 
photon emitted, as measured in the frame of the emit- 
ting matter. (1 + 0cos(6)) denotes an average weighted 
by the received intensity. We shall use the highly rela- 
tivistic limit — > 1 throughout this work, which leads to 
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the result (1 + /3cos(0)) = 4/3 applied in the final lines of 
equation |22| (Wijers & Galama 1998). 

We estimate the peak flux density following equations 19 
and 25 of Wijers & Galama (1998). The basic equation is 



r x N e x 



V3e 3 B 1 + z 

x 



p „, i 
m e c 



n^d 2 



(23) 



Here (j) p \/3e 3 B / (m e c 2 ) is the average comoving frame peak 
luminosity per unit frequency emitted by a single elec- 
tron. The details of the average over pitch angle and elec- 
tron energy are hidden in the factor <j> p , which is a func- 
tion of the electron energy distribution index p with range 
0.59 £ <j) p £ 0.66 for 2 < p < 3 (Wijers & Galama 1998). 
The factor T accounts for the Lorentz transformation of 
flux density (Wijers & Galama 1998). Distance and beam- 
ing effects enter through the factor 1/(£! 7 gP), where d is 
the luminosity distance to the burst, and f2 7 is the solid 
angle into which radiation is beamed. Finally, redshift af- 
fects the flux density by factor (l + z) (e.g. Weedman 
1986). Our goal is to express F v _ m ^ purely in terms of 
the dynamical variables we calculated in section |^. 
The comoving frame synchrotron cooling time is 



t. 



67T77l e C 



67rm e c 



u T ^ e B 2 <j T 7 e ,pcakS i 



(24) 



where o T = 0.665 x 10 _24 cm 2 is the Thompson cross- 
section (e.g., Rybicki & Lightman 1979). 

To obtain the magnetic field strength B, we need the 
volume of the ejecta cloud, which will have transverse ra- 
dius ~ r£ m + c s t co and thickness ~ c s t co , giving comoving 
volume V = ir{ctC, m + c s t co ) 2 (c s t co ). Under the approx- 
imation of negligible radiative losses, the internal energy 
is given by Ei tCO = Eq/T. The comoving frame magnetic 
field strength is thus 



B = 



1/2 



(25) 



The remaining pieces are trivial. fi 7 « n (£ m + 1/r) , 
and d and 1 + z are simply scale factors. 

We now have all the pieces of equations ^ and ^3] ex- 
pressed in terms of dynamical variables from section |[ 
This means that we can insert these formulae into our nu- 
merical integration code and calculate F v ^ m ^ and ^©, m as 
a function of t@ (or of /, T, or r). In order to determine 
a light curve at fixed observed frequency, we combine the 
broken power law spectral shape described above with the 
calculated frequency and flux density of the spectral peak 
to determine the approximate flux density at the observed 
frequency and time. 

The cooling break and self absorption break (cf. Sari, 
Piran, & Narayan 1998) are additional observed features 
in afterglow data. We do not treat either in detail here, 
but do we present a derivation of the cooling break be- 
havior for beamed gamma ray bursts elsewhere (Rhoads 
1999b). These results are summarized below. We have 
not yet treated the self-absorption break. Self-absorption 
is important primarily at low frequencies, where scintil- 
lation can hamper light curve slope measurements. For a 
treatment of this regime in beamed bursts, see Sari, Piran, 
& Halpern (1999). 



Finally, we consider the evolution of the apparent angu- 
lar size 9. In the spherical case or the power-law regime 

for a beamed burst, 9 — r/(Tdg) oc i^ 8 (where dg is the 
angular diameter distance to the burst). In the exponen- 
tial regime, 9 is determined by the physical transverse size 
of the ejecta cloud rather than the beaming angle, but 
the difference is not dramatic because the physical size in- 
creases at c s ~ c. The result is therefore 9 w c s t co /dg oc 

1 /2 

1/r oc ijg . There is also an intermediate regime, valid for 
the brief time when 1/r ^ Cm ^ {c s t co )/ (ct). In this case, 

9 oc £ m r oc t 1 ^ 4 . If the exponential regime did not happen 
at all, this behavior would continue for all T < l/£ m . 

3.2. Analytic Results: Optically thin case 

In the limiting cases where one of the terms in the trans- 
verse size ( m ct+c s t co is dominant and the other negligible, 
we can derive analytic expressions for F um ^ and ^© jm as 
functions of observed time t® and the physical parameters 
of the fireball. We begin with the early time case, and 
show that its light curve is observationally indistinguish- 
able from that of an isotropic burst. 

3.2.1. Power Law Regime 

We first determine the comoving magnetic field in this 
regime by inserting (, m ct 3> c s t co into equation |2^ to ob- 
tain 



B = 



3 3/8 1/2 ? B 



p(l + z)\ 3/8 



Xm7 V t 

Inserting this result into equation |2^, we find 

2 l/4 5 l/2 



(26) 



3 3 /87rV 



l/2 X p£e£ 



2*1/2 



m. 



1/8 V2 



P \ 1 / 8 /„ \ 3/8 

Eo\ (p_\ r3(1 



irC 2 

5V2 

2 7 /2 



>-5/8 



(27) 



x £ 2 £ 1/2 



m,, 
m, 



2 !/2 

m e c z c s 



— Y /2 r 3/2 n +z) 1/2 



(28) 



where we have used equation [IT] to eliminate T in the last 
line. 

For the cooling break, we obtain z/© iC ooi oc t ffl (Rhoads 
1999b; Sari et al 1999). 

Turning our attention to i^m,©, we first need the num- 
ber N e of radiating electrons in terms of £® . For the power 
law regime, this becomes 



N K = 



2 3 / 2 Hip 
3V4 M , 



1/4 



Eat 



Ot© 



cl + z) 



3/4 



(29) 



Combining this with equations |23|, |26|, and the appropriate 
limiting form of f2 7 , we find 



c p 1 / 2 E Q 1 + z 




p, e m v m e c 



(30) 
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Note that this result is independent of t®. 

Apart from small differences in the numerical coeffi- 
cients, our results for the power law regime are essentially 
the same as the results that Waxman (1997a, b) and Wi- 
jers and Galama (1998) obtained for isotropic bursts. Dif- 
ferences between our results and Waxman's are primarily 
because we have adopted the more precise treatment of 
the synchrotron peak frequency presented by Wijers and 
Galama (1998), while differences between our results and 
those of Wijers and Galama stem from a slightly different 
way of calculating the comoving frame magnetic field. 



with equations g3| and 32 to obtain 
Z l ' 2 h 2 



TP. 



1/2. 



E 3/2 C 2 

^0 w 



m e m p ^ e c 2 c 2 



x (c s W)- 3/2 f(^ _I 



1 + 2 
d 2 



(38) 
(39) 



Substituting the initial conditions for the exponential 
regime, this becomes 



3.2.2. Exponential Regime 

When c s t co Cmd, we are in the regime where T, tm , 
etc. all behave exponentially with radius (section 2.2.2). 
We first rewrite the scalings from equation pG as 

WW = exp{2(r-r fc )/r r } 

f/h = exp{2(r - r b )/r r } = t m /t m , b 
T/T b = exp{-(r-r b )/r r } = (W*©,i>)~ 1/2 
tco/t co ,b = exp{(r - n)/r r } = (t©/t e . b ) +1/2 • 

(31) 

We next determine the comoving frame magnetic field 
in the appropriate limit: 



B = 



1/2 



p-1/2 -3/2 
1 b b co,b 



-1/2 



(32) 



Combining this result with equation g2| allows us to deter- 
mine the peak frequency as 



1 2 3 / 2 



1 + Z 7T 

p l/2p5/2 r v_ 3 /2 
X ft 1 b \ c stco,b) 



V 


2 


e 






m e c 


t 


B_ 


r 


k 


,b 





(33) 



Substituting for the exponential regime initial conditions 
from equations then yields 



1 



2 ll/2 



T t2tl/2 I ™P 



m e \cJ Cm \t®,b 

ol/6cll/6 / m 

(l + z) 6 5 X f 2 £V2 ( ?h 



(34) 



(I) 



3/2 E, 



2/3 



1/6 10/3 



(35) 



P 



The observed cooling break frequency ceases to evolve in 
this regime: ^©, CO oi oc t^ (Rhoads 1999b; Sari et al 1999). 

Turning now to the amplitude of the spectral peak, we 
combine 



25 EoC 2 



8 c 2 s [i e m v t® jb 



and 



7r r b (WW) 



(36) 



(37) 



F„ 



(?) 



Bis f i/>_fl 
125 rm p m e ft v c 

-g yp ft® y 1 1 + . 



3/2 



(40) 



3 l/3 5 7/6 



1/2. 



/r \ 1/2 

(7) (41) 



2 7/2 7r 5/6^fl c 3 mpTOeA i £ V C 

£ 4/ y /6 (i + ^) V 



5/3 



d 2 



At tm = tmh , equations |34| and [lO] differ from equa- 
tions |2^ and 30 by factors of order unity. This difference 
is not worrying since our analytic approximations are not 
expected to be particularly accurate in the transition be- 
tween the two limiting cases. Numerical correction factors 
to the coefficients of equations U and ^ can be derived 
from n umeric al integrations. Such factors are presented in 
section 3.2. 5| below. 



3.2.3. TV Dinner Equations 

We now pause a moment to consolidate our results so far 
and express the key equations in terms of fiducial param- 
eter values 2 . We begin with equations ^8] and |3^. These 
become 



i/ ffl , m = 9.6x 10 12 x {l + z) 1 ' 2 



V0.525/ V 0.1 / V 0.1 



e/4 



ffo/10 53 erg \ 1/2 (t®_\ 



c/V3 
1/2 

-3/2 



-1/2 



(42) 



V da y/ 



and 



11 x (l + z) 

1/2 



0.1 



p 



0.637 Vc/V3 
1.3\ /£ /10 53 erg 

Ve) \ C m /4 
1/2 / d 



Hz 



-1/2 



(43) 



mjy 



2 We call these "TV dinner equations" because numerical values for 
further preparation. 



10- 24 g/cmV V 4 - 82G P C , 

physical constants have been inserted, so they are ready to use without 
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The observed time corresponding to the transition between 
the power law and exponential regimes is 



i ffi>6 = 12.1 x(l- 
P 



x 



10- 24 g/ cm 3 



-1/3 



~ 8/3 ^ /10 53 erg^ 1/3 



C/4 



8/3 



or I days 



(44) 



Thereafter, the frequency and flux density at the spectral 
peak are characterized by equations ^8] and ^0|. Numerical 
integrations show that modest correction factors e v w 0.74 
and e F w 0.7 should be applied to these two equations at 
late times to compensate for ap proximations in the ini- 
tial conditions (see section 3.2.4| below). These have been 
incorporated in the following three equations. 

The observed frequency of the spectral peak at the time 
of the break is 



v ®,m.b 



i : x io 11 

l + z 

2 



( — 

V0.7 



0.1 



o.i 



74 0.525 
P 



c/V3 



7/2 



(45) 



10- 24 g/ cm 3 
The subsequent evolution is given by 



c 

o.i 



^ Hz 



,6 X 



and 



TP, 



= 0.41 x -4: 



Me 

x (l + z) 



0.7 0.63 

V10 53 erg/4 
d 



3/2 



(46) 



1/2 



1/2 



4.82 Gpc 



IO" 24 g/ cm 3 
l 



A- ) n,h . 



Finally, for completeness, we include our result for 
^e.cooi from Rhoads 1999b: 



^ffi.cool — 



5.89 x 10 13 (*©/*©,&) 1/2 + 1.34 xlO 1 

17/6 / c \ -3/2 



Hz 



1 



l + z 
f p ■ cm 



c/V3 



0.1 



(48) 



: - 5/6 / J B /10 53 erg\- 2/ VC n 



V.10 



-24 , 



0.1 



-4/3 



Note that this equation already interpolates over the break 
time t©,f>; the int erpola tion was derived in the fashion sug- 
gested in section 3.2.5 below. 



3.2.4. Putting the Pieces Together 

An accurate description of the behavior in the transi- 
tion between the power law and exponential regimes can 
be obtained numerically. We first note that there is a sin- 
gle characteristic observed time i®^ given by equation |l8] 
and flux level i*i, )m ,e,& = Fv,m,®{t® < *©,b) given by equa- 
tion If we use these as our basic time and flux units, 



and denote the observed time and peak flux scaled to these 
units as f© and F„ imj ©, there is a unique Fi, !mi ©(t©) rela- 
tion. This is plotted in figure |^. 

Similarly, we can define the characteristic frequency 
v ®,m,b in the problem to be given by equation 28 eval- 
uated at t^fi, and v§~m to be the frequency scaled by 
this value. Then we can again obtain a unique relation 
is®^n(t®), which is shown in figure ||. 

At late times, the numerical integrations yield a flux 
density that is a factor e F 0.7 smaller than in equa- 
tion |40j, and a frequency of peak emission that is a factor 
e„ ~ 0.74 smaller than in equation |34|. This is presumably 
due to the approximate initial conditions used for the ex- 
ponential regime evolution. These initial conditions are 
obtained by applying an asymptotic approximation out- 
side its range of validity, and it should not be surprising 
if this procedure introduces some error. We suggest below 
that this error may be corrected empirically. 

To obtain predictions for a given set of model param- 
eters from these dimensionless curves, we need only (1) 
determine numerically the values of i©,&, -F^m,©,^ and 
^©,m,&; and (2) determine the time interval over which our 
assumption I/Tq <J / < Tq remains valid. The early be- 
havior, before the ejecta accrete a dynamically important 
amount of ambient medium (i.e., / < I/Tq) is unlikely to 
be observed at long wavelengths, since it is over within a 
fraction of a second for reasonable burst parameters. We 
therefore consider only the end condition, f ^ Tq. This 
happens at i© = i©,/, where 



~ IVffi,b/ fb 



25 



1 



Si) 



(49) 



At later times, our assumptions that T J> 2 and /3 « 1 
break down, and the behavior of the fireball changes again. 
Such changes may be relevant to the radio behavior of 
gamma ray burst afterglows, but we will not consider them 
here. 



3.2.5. Empirical Interpolations 

To obtain a readily calculated burst behavior around 
time t©,&, we can interpolate between the asymptotic be- 
haviors for earlier and later times. We do this first for 
Fv,m,<s and then for i^©, m . We use interpolants of the 
form g — (<7^ K + e$^ re ) _1 / K , where g\ and eg2 represent 
limiting behaviors of an arbitrary function g for early and 
late times. The exponent k determines the smoothness of 
the transition between the limiting behaviors. The scalar 
e is introduced so that the numerically derived correction 
factors to the late-time asymptotic results can be applied. 

For F„ )m ,©, the asymptotic behaviors are F„ )m ,© con- 
stant and F^ra,© ~ i© . We work with the scaled quan- 
tities defined in section fo.2.4 , so that the break between 
the two asymptotic behaviors is expected for log(i©) ~ 0. 
We set g\ — F u _ m ^^- We use equation ^ for g 2 , and set 
the correction factor e = 0.7. Finally, we choose n = 0.4. 
The resulting interpolation is plotted atop the numerical 
integration results in figure |^. 

The asymptotic behaviors of z^©, m are oc £© 3 ^ 2 and 
oc t^ 2 . In this case, we have taken g\ from equation p8[ 
For g 2 , we take equation 53, and set e = 0.74. Here we 
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find k = 5/6 works well. This interpolation is shown in 
figure HI 

3.2.6. Light Curves: Optically thin case 

The afterglow light curve at fixed observed frequency 
is obtained by combining the predicted behavior of f® lTTl 
and -F^m.® with the spectrum for a truncated power law 
electron e nerg y distribution. We use the analytic results 



Y 2 C, 2 m once r < i/c 



of section 3.2. Then we find four generic behaviors, de- 
pending on whether the frequency is above or below v^^ m 
and whether the time is earlier or later than These 
are 



F v ,® oc 



1/2 

3(P- 
-1/3 



l)/4 



*s 


< k 


3,6 ; 




< "a 


,m(k 


5) 


*a 


< k 


5,6 ; 




> "a 


,m[k 


5) 


*e 


> k 


3,6 ; 




< "a 


,m{k 


5) 


*e 


> k 


5,6 ; 


fa 


> 


,m(k 


5) 



(50) 

Here p is the electron energy spectrum slope and a = 
(p — l)/2 is the high frequency spectral slope, as usual. 
Note particularly how steep the light curve becomes for 
t® > t 9tb and i/ @ > u @tm (t @ ). 

Three representative light curves are shown in figure |4[ 
These have been derived by combining the empirically in- 
terpolated v^/m, and F vm j$ curves with the broken power 
law spectrum. Note that the rollover at the beaming tran- 
sition (log(f©) ~ 0) is rather slow, so that observed be- 
havior will be intermediate between the asymptotic power 
laws of equations p0| for a considerable time. This slow 
rollover is in part aue to the compound nature of the 
break. The light curve decay begins to accelerate as soon 
as we can "see" the edge of the jet, when T < 1/Cm- The 
additional steepening when dynamical effects of beaming 
become important occurs slightly later, when r < Tf, ~ 
0.23/Cm (cf. equation [l7[) (cf. Panaitescu & Meszarosl999 
for additional discussion of this point). 

Equation |5C| assumes ^©, a bs < ^® < ^©,cooi throughout. 
(Here ^® ja bs is the self-absorption frequency, measured in 
the observer's frame.) If we now include the cooling break, 
we obtain the additional light curve behaviors (derived in 
Rhoads 1999b) 



F u 



oc 



,l/2-3p/4 

k 

ten to 



< k 
> k 



,cool(^ffi) 
,cool(i©) 



(51) 

where we have also assumed that ^© >m < v®,coo\- 

These 

behaviors are not shown in figure^, but were used to fit the 
light curve of GRB 970508 with beamed afterglow models 
(Rhoads 1999b). 

3.2.7. Light Curves: Optically thin case without sideways 

expansion 

It remains possible to constrain gamma ray burst beam- 
ing by looking for light curve breaks even in the case where 
lateral expansion of the evolving burst remnant is unim- 
portant. This corresponds to dropping our fourth assump- 
tion from section ||. 

In this case, th e dyn amical evolution follows power law 
behavior (section 2.2.1) throughout, but the emergent ra- 
diation is diluted relative to the spherical case by a factor 



V 2 oc t, 



-3/4 



For an adiabatic evolution, we have 
The power law exponents for the light curve 



in this regime become —1/4 for ^© ja b s < v < v®, m , ~3p/4 
for z/© )Tn < v < ^e,cooi, and -1/4 - 3p/4 for ^©, CO oi < v. 

The most plausible mechanism for quenching lateral ex- 
pansion is a strongly radiative remnant, so it is more 
relevant to examine this regime with radiative dynam- 
ics. Then T 2 oc £© 6 ^ 7 , and all the light curve expo- 
nents decrease by an additional 3/28, becoming —5/14 for 
^®,abs < v < v® >m , -3/28 - 3p/4 for v e . m < v < v 9<coo i, 
and —5/14 — 3p/4 for ^© jC ooi < v. 

While t hese s lope changes are less dramatic than those 
in section 3.2.6f they would be strong enough to detect in 
afterglow light curves with reasonably large time coverage 
and good photometric accuracy. 

3.3. Optically Thick Case 

We now consider briefly the electron energy distribu- 
tion model of Paczyhski & Rhoads (1993). This model 
differs from that of the preceding sections in a few ways. 
First, the electron power law index was fixed at p = 2 
to avoid strong divergences in the total energy density in 
electrons. Second, the minimum electron energy £ min was 
taken to be sufficiently small that emission from electrons 
with £ = £ m i n was always in the optically thick regime. 
Under these circumstances, there is a single break in the 
electron energy spectrum at the frequency corresponding 
to optical depth r = 0.35 (cf. Pacholczyk 1970), with 
spectral slope v^l 2 below the break and v~ x l 2 above. The 
magnetic field behavior is the same in this model and more 
recent ones. 

Combining this electron behavior with the power law 
regime dynamical model reproduces the scalings from 
Paczyhski & Rhoads (1993), namely 



and 



cl/3 1/3.-2/3 
i/ e , m oc E ' p 1 t @ 



7/8 1/8, 5/8 ^ p 13/12 1/3,-5/12 



p _ oc E 1 n ' v 1 oc E n't 
r v,®,m oc rjQ p u ®.m ^ ^0 r 1 



(52) 



(53) 



If we use instead the exponential regime dynamical 
model, we find 



&.»>! 



OC t 



-1 



and 



F, 



v,®,m 



OC t 



-3/2 



(54) 



Readers interested in the precise numerical coefficients 
for these relations are referred to Paczyhski & Rhoads 
(1993) for the spherical case. For the beamed case, nu- 
merical results may be found by applying the Paczyhski & 
Rhoads (1993) results at the transition between the power 
law and exponential regimes, and continuing the evolution 
using equations |54|. 

The light curve for this electron model then becomes 



F v ,® oc 



5/4 
b ® 




< k 


3,6 ; 


^a 


< "a 


,m 


,-3/4 
b ® 


*a 


< k 


3,6 ; 




> v& 


,rn (t 




*G 


> k 


5,6 ; 




< fa 


,rn (t 


L ® 


*<E 


> k 


3,6 ; 


v& 


> fa 





(55) 



The behavior here at frequencies f© > v®,m{t®) is the 
same as in equations |5^ with p = 2. However, ^© im has a 
different meaning in the two models. 
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From these results, we see that the substantial changes 
in the observable behavior of a beamed burst are not de- 
pendent on the precise nature of the electron energy dis- 
tribution. 



4. DISCUSSION 

We now put mathematics aside to recapitulate our re- 
sults and to discuss their implications for the interpreta- 
tion of afterglow observations. 

We have shown that the dynamics of a gamma ray burst 
remnant change qualitatively when the remnant's Lorentz 
factor T drops below the reciprocal opening angle 1 / £ m of 
the ejecta. Before this time, the Lorentz factor behaves 
as a power law in radius. Afterwards, the Lorentz factor 
decays exponentially with radius. The change occurs be- 
cause lateral expansion of the ejecta cloud increases the 
rate at which additional material is accreted. Such lateral 
expansion is prohibited by symmetry in the spherical case. 

When the remnant enters this "exponential regime," the 
relation between the observed spectrum and the observed 
light curve changes. Inferences about the electron energy 
spectrum in afterglows come from the light curve decay 
rate and spectral slope. The general agreement between 
the two methods has been taken as a confirmation of the 
(spherically symmetric) fireball model (Wijers, Rees, & 
Meszaros 1997; Waxman 1997a). 

The light curve decline at frequencies above the spec- 
tral peak becomes very steep (i® p , where p is the index of 
the electron energy spectrum) once the burst dynamics en- 
ter the exponential regime. Reconciling this relation with 
the observed decays (—1 £ d\og(F u ^)/d\og(tQ) £ —1.5) 
would require an extremely flat electron energy spectrum, 
and consequently a very blue spectral energy distribution. 
This was not seen in early observed spectral energy distri- 
butions (see Wijers et al 1997 for GRB 970228; Reichart 
1997 and Sokolov et al 1997 for GRB 970508; and Re- 
ichart 1998 for GRB 971214). We infer that GRBs 970228, 
970508, and 971214 were probably not in the exponen- 
tial regime during their observed optical afterglows. GRB 
971227 provides a possible, though dubious, counterexam- 
ple. There is one image suggesting a counterpart (magni- 
tude R 19.5) on December 27.91 (Castro- Tirado et al 
1997). Later images show no corresponding source, requir- 
ing a decay at least as fast as i® 2 ' 5 (Djorgovski et al 1998). 

This is consistent with a i® p decay for typical values of p. 
However, this explanation remains speculative, since there 
is no second image confirming the proposed counterpart. 

Subsequent afterglows have provided a more hopeful pic- 
ture for practical application of beaming models. In par- 
ticular, GRB 990123 shows a break that is quite possibly 
due to beaming (e.g., Castro-Tirado et al 1999; Kulkarni 
et al 1999), and comparison of the spectral slope and decay 
slope for GRB 980519 gives better agreement for beamed 
than for spherical regime models (Sari et al 1999). 

In the case of GRB 970508, we can place a stringent 
limit on the beaming angle in the context of our model. 
The optical light curve extends to ~ 100 days after the 
burst and does not depart drastically from a single power 
law after day 2 (Pedersen et al 1998); thus, no transition 
to the exponential regime occurred during this time. As 
already noted, the spectral slope and light curve law decay 



rate are in fair agreement for the spherical case, and poor 
agreement for the beamed case. The radio light curve 
furnishes the last critical ingredient. Goodman (1997) 
pointed out that diffractive scintillation by the Galactic 
interstellar medium is expected in early time radio data, 
and that this scintillation will stop when the afterglow 
passes a critical angular size. By comparing this charac- 
teristic size with the time required for the scintillations to 
die out, one can measure the burst's expansion rate. This 
test has been applied (Frail et al 1997; Waxman, Kulka- 
rni, & Frail 1998) and shows that r <; 2 at f ® ~ 14 days. 
Thus, no power-law break to faster decline is observed at 
r ;> 2, and we infer that GRB 970508 was effectively un- 
beamed (£ m > 1/2). This rough derivation is borne out 
by detailed fitting of beamed afterglow models to the GRB 
970508 light curve, which yields the same beaming limit 
C,„ £1/2 radian (Rhoads 1999b). 

This conclusion, combined with the GRB 970508 red- 
shift limit z > 0.835 (Metzger et al 1997), immediately 
implies a minimum energy for the burst. This burst was 
detected as BATSE trigger 6225, and the total BATSE 
fiuence was (3.1 ±0.2) x 10 -6 erg/ cm -2 over the range 20- 
1000 keV (Kouveliotou et al 1997). The gamma ray emis- 
sion alone therefore implies Eq £ 4.7 x 10 51 (f2/47r) erg £ 
3 x 10 50 erg. Here we have based the luminosity distance 
on an Q = 0.2, A = 0, H = 70km/s/Mpc cosmology, 
and applied the beaming angle limit £ m ;> 0.5 radian in 
the second inequality. This conclusion is of course model- 
dependent and might change if our assumptions about the 
blast wave physics or beaming geometry are badly wrong. 
We will discuss possible ways to reduce the energy require- 
ments of GRB 970508 wh ile r etaining consistency with the 
afterglow data in section 4.1 below. 



If the beaming angle £ m is substantially variable from 
burst to burst, it is possible that some bursts enter the 
rapid decay phase before the spectral peak passes through 
optical wavelengths. Present data suggests that this is 
indeed the case; GRB 980519 is best fit by assuming ex- 
ponential regime behavior (Sari et al 1999), while GRB 
990123 appears to be a transition case with a break ob- 
served in the optical light curve (e.g. Castro-Tirado et al 
1999, Kulkarni et al 1999, Sari et al 1999). The resulting 
rapid decay could then explain some of the optical non- 
detections of well studied GRBs such as 970828 (Groot et 
al 1997). Alternatively, for characteristic beaming angles 
1 C™ £ 0.1, we would expect beaming to become dy- 
namically important between the time of peak optical and 
radio afterglow. This would then help explain the paucity 
of radio afterglows, which unlike optical afterglows cannot 
be hidden by dust in the burster's environment. There is 
some evidence that the radio emission involves a different 
process, or at least a different electron population, from 
the optical and X-ray afterglows: The peak flux density in 
GRB 970508 did not follow a single power law with wave- 
length as it ought to under the simplest fireball models 
(Katz & Piran 1997b). 

The transition in light curve behavior at T ~ 1/Cm 
is also important for "blind" afterglow searches. Such 
searches would look for afterglows not associated with de- 
tected gamma ray emission. A much higher event rate 
for afterglows than for bursts is a natural consequence of 
beamed fireball models, since the afterglow emission peaks 
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at lower bulk Lorentz factors than the gamma ray emis- 
sion does. Comparison of event rates at different wave- 
lengths can therefore constrain the ratio of beaming an- 
gles at those wavelengths (Rhoads 1997a). However, we 
will only see the afterglow if either (a) we are within angle 
Q m of the burst's symmetry axis, and therefore could also 
see the gamma ray burst, or (b) the Lorentz factor has 
decayed to V < 1/Cm and the afterglow light curve has 
entered its steep decay phase. We have already argued 
that GRB 970228, GRB 970508, and GRB 971214 were 
not in this steep decay phase based on the comparison of 
light curves and spectral slopes. It follows that if blind 
afterglow searches find a population of afterglows not as- 
sociated with observed gamma rays, those afterglows will 
exhibit a steeper light curve decay than did the 1997 af- 
terglows. The efficiency for detecting such rapidly fading 
"orphan" afterglows will be substantially lower than the 
efficiency estimated from direct comparison with spherical- 
regime afterglows. 

Other models of beamed gamma ray bursts are possible. 
In particular, we have assumed a "hard-edged" jet, where 
the mass and energy emitted per unit solid angle are con- 
stant at small angles and drop to zero as a step function 
at large angles. Profiles in which these quantities decrease 
smoothly to zero may be more realistic. Whether these 
differ importantly from the model presented here depends 
on whether most of the energy is emitted into a central 
core whose properties vary slowly across the core. Lay- 
ered jet models in which most of the kinetic energy from 
the fireball is carried by material with a low Lorentz factor 
can have substantially different afterglow light curves from 
either the spherically symmetric case or the hard-edged jet 
case. This is because the afterglow emission can be dom- 
inated by outer layers where the initial Lorentz factor is 
high enough to yield optical emission during ejecta deceler- 
ation, but insufficient to yield gamma rays. The afterglow 
is thereby effectively decoupled from the gamma ray emis- 
sion, and it becomes harder to predict one from the other. 
Such models have been explored by several groups (e.g., 
Meszarosfe Rees 1997b; Meszaros, Rees, & Wijers 1997; 
Paczyhski 1997). A similar decoupling of the gamma-ray 
and afterglow properties can be produced in the spherical 
case by allowing inner shells of lower Lorentz factor and 
larger total mass and energy to follow the initial high-r 
ejecta (Rees & Meszaros 1998). 

It is possible to approximate the afterglow from a lay- 
ered jet by a superposition of hard-edged jets. For this 
to be reasonably accurate, the outer layers should have 
Lorentz factors substantially below those of the inner lay- 
ers, and opening angles and energies substantially above 
those of the inner layers. 



4.1. Energy Requirements for GRB 970508 

We now consider how our model will change if we vary 
some of the basic assumptions. Our primary concern is to 
determine whether the minimum energy required to power 
the GRB 970508 afterglow can be reduced substantially 
below the requirements derived from a spherical adiabatic 
fireball model expanding into a homogeneous medium. We 
will therefore sometimes err on the side of extreme model 



assumptions chosen to minimize the energy needs. In or- 
der to declare a model consistent with the data, we require 
that either (1) there be no break in the light curve or spec- 
trum around T ~ 1/Cm, or (2) the break occurs early (be- 
fore t® ~ 2 days) and the late time light curve shows a slow 
decline even for spectral slopes as red as those observed. 

The first requirement is physically implausible. Even 
in the absence of the dynamical effects reported above, so 
long as the afterglow is from relativistically moving and 
decelerating material, its flux will scale with an extra fac- 
tor T 2 once F ^ 1/Cm- Since T decreases with time, a 
break is generally expected, though perhaps it could be 
avoided with sufficient fine-tuning of the model. 

The second possibility is more interesting. It requires 
us to construct a model where factors besides beaming 
contribute relatively little to the decay of F v ^ m ^ with t®, 
or where the observed spectrum does not directly tell us 
about the electron energy distribution. A burst expand- 
ing into a cavity (such that p increases with r) might 
give a slow decay, while a sufficiently large dust column 
density would give a red spectrum despite a flat electron 
energy distribution (cf. Rcichart 1997, 1998). However, 
both would require some degree of fine tuning. The dust- 
reddened spectra would deviate measurably from pure 
power-laws given good enough data, but the present data 
are probably equally consistent with both pure and red- 
dened power law spectra. Certainly such reddened beam- 
ing models would imply little correlation between observed 
spectral slopes and light curve decays, since the dust col- 
umn density could vary wildly from burst to burst. This 
hypothesis is somewhat ad hoc, but is consistent with 
present data and is supported by other circumstantial evi- 
dence linking GRBs to dust and star forming regions (e.g. 
Groot et al 1997; Paczyhski 1998). At present, then, it ap- 
pears the most viable way of reconciling beamed fireball 
scenarios with the 1997 afterglow data. 

We now discuss a few variations of fireball models in 
greater detail. 

4.1.1. Radiative case 

We first consider the behavior of a radiative regime fire- 
ball. In this regime, the internal energy of the fireball is 
low, since it is converted to photons and radiated away. 
The largest implications for beaming are when the inter- 
nal energy density is so low that c s <C c. In this case, 
the lateral expansion that leads to the exponential regime 
of burst remnant evolution in the adiabatic case is unim- 
portantly small. We assume this low sound speed through 
much of the following discussion. 

We assume that energy in magnetic fields and protons is 
transferred to electrons in the burst remnant on a remnant 
crossing time (<~ t co ). The electrons are assumed to main- 
tain a power law energy spectrum, with a large £ ma x whose 
precise value is determined by the requirement that the 
burst radiate its internal energy efficiently. Under these 
circumstances, the Lorentz factor scales as F oc r~ 3 and 
the comoving frame internal energy Ei nt of the remnant 
follows the evolution 

dE int 2 2 2 E int alt co npc 2 C^(Tr 3 ) 4E lnt 

— ~r pc7 rc m r ~—^r — • 

(56) 
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This admits a solution of the form Ei nt ~ (nC^pc 2 ^^) — 
C2^~ 4 )/4 where C2 is a constant of integration. At late 
times, we throw away the cir~ term, which becomes neg- 
ligible. The result then becomes 



A = 0, and H = 70km/s/Mpc (cf. equation M). We 
compare this to the optical fluence of the burst, which we 
estimate as 



E int « (Tr 3 )Trt 2 mP c 2 /4 



(57) 



Qopt — 



which is constant since T 



„-3 



in this regime. If the 



o Jo 



Fvmdvdtq 



(59) 



sound speed becomes negligibly small at some point in 
the burst remnant evolution, then the volume of the shell 
scales as V oc r 2 thereafter. The magnetic field then scales 
as B cx r , based on constant Ei nt . The observed peak 
frequency scales as ^©, m cx T 3 B cx r~ 10 cx t^ 10 ^ 7 . 

The power radiated is simply ~ T x (r — l)c 2 x ■nQ n r 1 cp 
in the observer's frame, but this is dominated by emis- 
sion from electrons at £ ma x, which we have not calcu- 
lated. The peak in F v q must be estimated as before. We 
have total comoving frame energy ~ £, e Ei n t in electrons at 
£ ~ £ m inj which is radiated over the comoving cooling time 
t s ~ l/(TB 2 ) ~ V/(TE int ) ~ r 5 . In the observer frame, 
this gives total power output ~ T 2 Ei nt /t s , accounting for 
factors of T from the Lorentz boost to photon energies and 
for the transformation between t co and t& . The frequency 
range containing this power scales as Az/© cx ^©, m . So, 
in the spherical case, i^m,© cx T 2 Ei nt /(t s v^, tm ) cx r^ 1 cx 

— 1/7 

to, .If we now allow for beaming, we introduce another 



Inserting our broken power law spectrum and the depen- 
dence of ^©, m on i© yields 



! opt 



p-7/3 



1/3 



factor of fi 7 1 



T 2 /ir and obtain F v<mt ® cx r 7 act 



( 6 °) 

Here v\ is an arbitrary reference frequency, and i®,m(fi) 
is defined as the moment when ^® lTO = v\. Setting 
v\ = 6 x 10 14 Hz (corresponding to wavelength 0.5/zm), 
^,m,® = 30/uJy, t©, m (^i) = 2 days, and p » 2.85 (cor- 
responding to a ' light curve) yields Q op t = 3.0 x 
10 _7 (^ ma2; /^i) 1 / 3 erg/ cm 2 . (A similar calculation using 
p = 2.2 and accounting for the additional break at the cool- 
ing frequency yields a similar fluence, 3.8 x 10~ 7 erg/ cm 2 , 
for v ma x — v \ — 6 x 10 14 Hz. This is the value used 
in Rhoads 1999b.) Taking luminosity distance 4.82 Gpc 
and considering only optical and longer wavelength after- 
glow, the smaller optical fluence estimate implies E 



Finally, putting in the spectral shape for fixed > 4 - 5 x 10 50 (fi/47r) erg £ 2.8 x 10 9 erg. We have applied 



we find that F v , 



cx i e ^ 5p+2 ^ 7 . Thus, this radiative 
regime model yields scalings fairly similar to our canonical 
adiabatic model. In particular, the late time light curve 
again shows a steep decline. While the assumptions made 
here may not be fully self-consistent, allowing c s ~ c would 
likely further steepen this decline. 

This result suggests that the GRB 970508 data can- 
not easily be reconciled with a beamed radiative afterglow 
model. 

4.1.2. Beamed burst, isotropic afterglow 

Gamma ray bursters may give rise to both fast and slow 
ejecta, where "slow" here means too slow to cause gamma 
ray emission. In this case, the optical and 7-ray properties 
of the event may be effectively decoupled if the slow wind 
contains most of the energy. 

Suppose the gamma ray burst is caused by a small 
amount of extremely relativistic ejecta, while the after- 
glow is caused primarily by a greater mass of material with 
low To. The afterglow light curve places almost no direct 
constraint on the isotropy of the first (high-Fo) material. 
However, the low-r material must be reasonably isotropic 
to avoid a visible break in the light curve at late time. To 
explain a peak optical flux of ~ 30 pJy in the optical, we 
need total energy 



E a = 1.5 x 10 50 erg x 

/ 3\ -1/2 

/ p ■ crrr x 



0.63 p e 



1-24 , 



1.835 
1 + z 



1.3 



1/2 

(58) 



4.82 Gpc 



(see equations pG and |43j) , where we have assumed isotropy 
and where 4.82Gpc is the luminosity distance correspond- 
ing to z = 0.835 for cosmological parameters 17 = 0.2, 



our beaming limit, ( m J> 0.5 radian, to derive the lower 
limit here. If we take v max corresponding to the soft 
X-ray afterglow, the energy rises by another factor of 
~ 10. These fluence-based energy needs are dangerously 
close to exceeding the energy requirements from equa- 
tion [5^. Since the latter equation is based on an energy- 
conserving model, this comparison shows that £ B must 
be substantially below 1 and/or the density substantially 
below 10 _24 g/cm 3 if the model is to be self-consistent. 
Otherwise, the total energy radiated is comparable to or 
greater than the total energy available. Reassuringly, the 
density and magnetic energy fraction found by Wijers & 
Galama (1998) are roughly consistent with these require- 
ments. This consistency check could be refined by replac- 
ing equation BCfl with a more detailed fluence calculation. 



4.1.3. Layered jet models 

In this class of models, considered (for example) 
by Meszaros, Rees, & Wijers (1997) and Panaitescu, 
Meszaros, & Rees (1998), the material dominating the 
emission changes continuously through a range of initial 
Lorentz factors. We can approximate such models as a 
superposition of many "hard-edged" jet models. We have 
tried developing such models while minimizing the energy 
~i/2requirements. To do this, we build a sequence of adiabatic 
hard-edged jets, enforcing either the condition ^© jTrl = v§ 
or the condition t^ t b = t® throughout the afterglow evo- 
lution, and then adjusting the input energy requirement 
to match the observed light curve. (Here denotes the 
fixed frequency at which our data was taken.) A prelimi- 
nary exploration of such models has not yielded any drastic 
reduction in energy requirements. A more thorough study 
may be needed to make this conclusion firm. 
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4.1.4. Inhomogeneous environments 

The predictions of fireball models change somewhat if 
the ambient medium is not uniform. To date, investiga- 
tions of variable density environments have concentrated 
on density laws p oc r~ 9 (Vietri 1997b; Meszaros, Rees, 
& Wijers 1997; Panaitescu, Meszaros, & Rees 1998). The 
best motivated choices of g are g = and g = 2, which 
correspond to a uniform density medium and the den- 
sity profile expected from a constant speed wind from the 
burst progenitor expanding into a vacuum (or an ambient 
medium of much lower density) . When the ambient den- 
sity decreases with increasing distance from the burster, 
the general result for spherical symmetry is a faster decay 
of the afterglow flux (e.g. Panaitescu et al 1998) (though 
the duration of the afterglow should increase correspond- 
ingly). We therefore infer that a decreasing density profile 
will also steepen the light curve decline in the beamed 
case. This only exacerbates the disagreement between the 
observed slow afterglow decay and the model predictions 
for beamed bursts, given the observed spectral energy dis- 
tribution. 

It is worth asking how t he exp onential regime of burst 
remnant evolution (section 2.2.2] ) will change if the ambi- 
ent density is nonuniform. The exponential scale length 



oc p 



-1/3 



in the uniform density case. We therefore 



conjecture that a solution similar to the following may be 
possible: 



r 



so that 



exp 
exp 

r oc exp 



r T (Pi) x (P/Pi)~ 1/3 
r 

ry(pi) x (r/ri)s/ 3 

r \ (3-ff)/3' 
'V 



(61) 



(62) 



where r% and p\ are some fiducial radius and the corre- 
sponding density, and where r\, = (r r (pi) 3 r^ 9 ) ^ 3 9 \ 



4.1.5. Observational 



concerns 



Throughout this discussion, we have tacitly assumed 
that the spherical fireball does fit the GRB 970508 obser- 
vations well, so that difficulties with the various beaming 
models offer support to the spherical model. This is open 
to question. In defense of the spherical model, Reichart 
(1997) has studied the largest available afterglow data set 
from a single optical observatory (that of Sokolov et al 
1997) and finds that the afterglow of GRB 970508 is well 
fitted using a standard model with the addition of a mod- 
est amount of dust extinction at high redshift. By using a 
single data set, he minimizes many of the possible system- 
atic errors, such as inconsistent zero points for absolute 
photometry. 

On the other hand, if one examines the spectral slope 
from mixed data sets over larger wavelength intervals (op- 
tical - near infrared) and a larger time range, we find some 
worrying data points. In particular, the HST observations 
(Pian et al 1998) exhibit a spectral slope a w 1.5 ± 0.3 
based on quoted R (0.7pm) and H (1.6pm) band mag- 
nitudes from the STIS and NICMOS instruments. The 



observed slope in Sokolov et al's data is 1.1, which Re- 
ichart interprets as a reddened spectrum with intrinsic 
slope 0.8. The HST data point is thus in mild conflict 
with the Sokolov et al observation. The significance of 
this conflict is unclear, since the error on the HST data is 
dominated by calibration uncertainties. 

Likewise, if we compare the spectral slopes inferred from 
the Keck K s band data (Morris et al 1997) and near con- 
temporaneous optical data (Djorgovski et al 1997; Kele- 
man 1997), we find slopes of 0.24 ± 0.12 at t m = 4.35 day 
and 0.40 ± 0.10 at = 7.35 day. These are now substan- 
tially bluer than the value from Sokolov et al. The first of 
these may simply indicate that the K band flux has not 
yet passed its peak and entered the power law decay phase; 
scaling from the R band peak at i® « 1.9 day gives a K 
band peak at ~ 4.1 day. The second is harder to explain 
physically but easier observationally, because the R band 
flux estimate is based on an unfiltered observation that 
may have (for example) a substantial color term. 

The net effect of such outlying data points is illustrated 
by the light curve fits in Rhoads 1999b. These fits achieve 
X 2 per degree of freedom around 3.6 in fitting to a large 
compilation of R band data (Garcia et al 1998). It is not 
likely that any current model can do better without dis- 
carding either predictive power or outlying data points. 

In summary, the spherical model does fit the GRB 
970508 afterglow model better than the beamed model 
developed in this paper for any beaming angle £ m < 0.5 
radian. There are a few discrepancies in the measured 
spectral slopes. If these are real, they pose a challenge to 
standard fireball models. However, they could merely be 
indicative of calibration problems in inhomogeneous data. 
It is noteworthy that Sokolov et al (1997), who have the 
largest multiband data set from a single telescope, find 
no evidence for spectral slope evolution over the interval 
2 day < t© < 5 day. 



5. CONCLUSIONS 

We have shown that under a simple model of beamed 
gamma-ray bursts, the dynamical evolution of the burst 
remnant changes at late times. This change introduces 
a break in the light curve, which is potentially observ- 
able. The afterglows of GRB 970508, 970228, and 971214 
showed no convincing evidence for such breaks, and their 
combined spectral slopes and light curves are inconsistent 
with the predictions of this beamed model. This implies 
that beaming tests based on blind searches for afterglows 
must be prepared to identify transients whose properties 
differ appreciably from the properties of these approxi- 
mately spherical afterglows. Some more recent afterglows 
(GRB 990123; GRB 980519) better match beamed burst 
models, and may provide more suitable templates for these 
searches. 

Comparing our model with late time optical and ra- 
dio observations, we suggest that GRB 970508 was effec- 
tively unbeamed. This implies energy requirements that 
are near the canonical isotropic values for cosmological dis- 
tances, and are not greatly mitigated by strong beaming. 
No straightforward variation on our beamed fireball model 
seems likely to simultaneously explain the observed spec- 
tral slope and a pure power law light curve decaying at the 
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observed rate. We conjecture that strongly beamed fire- 
ball models cannot explain all observed gamma ray burst 
afterglows without substantially altering at least one ma- 
jor ingredient of the models. We therefore obtain the first 
lower bound on GRB energy requirements that does not in- 
volve assumptions about beaming: E ;> 3 x 10 49 erg. This 
limit will increase by an order of magnitude if the same 
material that gives rise to the optical afterglow causes ei- 
ther the X-ray afterglow or the gamma ray emission, and 
will rise further if the burst's energy is not converted to 
radiation with perfect efficiency. 

I wish to thank Ralph Wijers, Jonathan Katz, Tsvi Pi- 
ran, Eli Waxman, Daniel Reichart, Alexander Kopylov, 
David De Young, and Sangeeta Malhotra for helpful com- 
munications. I also wish to thank Infrared Processing and 
Analysis Center for hospitality during the course of this 
work. Finally, I wish to thank all those who worked to 
achieve accurate gamma ray burst position measurements 
and so opened the way for afterglow studies. This work 
was supported by a Kitt Peak Postdoctoral Fellowship. 
Kitt Peak National Observatory is part of the National 
Optical Astronomy Observatories, operated by the Asso- 
ciation of Universities for Research in Astronomy. 



14 



Submitted to The Astrophysical Journal, 1998; current revision 3/1999 



REFERENCES 



Bond, H E 1997, TAT.T Ckculai 6654 



Panmtr-srn , A 



Castro Tirado, A. J., Gorosabel, J., Grcincr, J., Zapatcro-Osorio, M. 



R., & Costa, E. lyy/, 1AU Circular 68UU 
Castro-Tirado, A. J., et al 1999, to appear in Science, 26 March 1999 
Costa, E. et al 1997, IAJJ Circular 657? , 
Dar, A. 1997, preprint ( |astro-ph/970923l| ) 

Djorgovski, S. C, Metzger, M. R., Kulkarni, S. R., Odewahn, S. C, 

Gal, R. R., and Pahre, M. A., Frail, D. A., Costa, E., and Feroci, 

M. 1997, IAU Circular 6660 
Djorgovski, S. G., Kulkarni, S. R., Ramaprakash, A. N., & Frail, D., 

on behalf of the Caltech GRB collaboration, 1998, GCN notice 

025 

Frail, Q A T^.lWrH B R Nir.rfrn T. PWnH M Xr TgyW CI 



B. 1 397, Nature, 389, 261 _ 

Garcia M. R. et al 1998, ApJ 500, L105 
Goodman, J. 1986, ApJ 308, L47 
Goodman, J. 1997, New Astronomy, 2, 449 
Groot, P. et al 1997, ApJ, 493, L27 
Katz, J. I. 1994, ApJ 422, 248 
Katz, J. I., & Piran, T. 1997a, ApJ 490, 772 
J I ,, U P , 



Katz, J. I., k Riran, T, 1997b, presentation at dth Huntsvillc GRB 

Syn posium, ^stro-ph/9712242 

Keleman, J. 1997, International Bulletin oj Variable Stars 4496. 
Kouvcliotou, C, Briggs, M. S., Preece, R., Fishman, G. J., Meegan, 

C. A., & Harmon, B. A., on behalf of the CGRO BATSE team, 

1997, IAU Circular 6660 . . 

Kulkarni, S. R., et al 1999, Nature, in press, |astro-ph/9902272 
Meegan, C. A., Pendleton, G. N., Briggs, M b., Kouvcliotou, C, 

Koshut, T. M., Lestrade, J. P., Paciesas, W. S., McCollough, M. 

L., Brainerd, J. J., Horack, J. M., Hakkila, J., Henze, W., Preece, 

R. D., Mallozi, R. S., & Fishman, G. J. 1996, ApJS 106, 65 
Morris, M., Mastrodemos, N., Zuckcrman, B., McCarthy, C, Beck- 

lin, E., Lowrance, P., Chary, R., and Barnbaum, C. lAU Circular 

6666 

J. 1997a, ApJ 476, 232 
J. 1997b, ApJ 482, L29 



Meszaros, P.. 
Meszaros, P.. 
Meszams, LL 



ph/ ^709273 
Metzger, ivi. R. , Djorgovski 



: Rees, M. 
: Rees, M. 
Rees, M. 



Wijers, R. A. M. J. 1997, astro- 



S G 
Frail, D. A., 



Kulkarni, S. R., Steidel, C. C, 
Costa, E., & Frontera, F. 1997, 



Adelberger, K. L 
Nature, 387, 879 

Pacholczyk, A. G. 1970, Radio Astrophysics, San Francisco: W. H 
Freeman 



ph/9801258 ) 
Fanaitescu, A." 



Meszaros, P., & Rees, M. J. 1998, preprint (astro- 



Meszaros, P., 1999 
1986, ApJ 308, L43 
, k, Rhoads, J. E 1993, 

1997, preprint (astro-ph/9706232|) 

1998, ApJ 494, T35 ^=1 

et al 1998, ApJ 496, ??? Qstro-ph/9710322 ) 



, submitted to ApJ 
ApJ 418,,L5 



Paczyhski, B. 
Paczyhski, B. 
Paczyhski, B. 
Paczyhski, B. 
Pedersen, H., 

Pian, E., Fruchter, A. S., Bergeron, L. E., Ihorsett, S. E.', Frontera 
F., Tavani, M., Costa, E., Feroci, M., Halpern, J., Lucas, R. A. 
Nicastro, L., Palazzi, E., Piro, L., Sparks, W., Castro-Tirado, A 

J., Gull, T., Hurley, K., & Pedersen, H. 1998, ApJ 492, L1Q3- 

-Efies, M 1 & Meszaros, P. 1998, ApJ, submitted (astro- 

ph/9712252) . 

Reichart, D. E . 1997, preprint ( astra=ph/aZl 2J nn 

Reichart, D. E. 1998, preprint ( astro-ph/9801139 

Rhoads, J. E. 1997a, ApJ 487, hr 

Rho ad i D. J. E, 1907b, 



presentation at 4th Huntsville GRB Symposium 



( |astro-ph/9712042D 
Rhoads, J. E. 1999b, to appear in A&AS (proceedings of the Rome , 

workshop oi l "Gamma Ray Bursts in the Afterglow Era"), astro- 

ph/9903400| (?) 

Kybicki, U. B., & Lightman, A. P. 1979, Radiative Processes in As- 
trophysics, New York: John Wiley Sz Sons 
Sari, R. 1997, ApJ 489, L37 

Sari, R., Piran, T., & Narayan, R. 1998, ApJ 497, L17 
Sari, R., Piran, T., & Halpern, J. P. 1999, astro-ph/99033 
Sokolov, V. V., Kopylov, A. I., Zharikov, S. V., Costa, E., Feroci, M., 
Nicastro, L., & Palazzi, E. 1997, presentation at 4th Huntsville 



GRB Symposium ( astro-ph/970909S ) 
Van Paradijs, J., et al lyy/, Nature 386, 686 
Vietri, M. 1997a, ApJ 478, L9. 
Vietri, M. 
Waxman, 
Waxman, 
Waxman, 
Wccdman. 



University Press 
Woods, E., & Loeb, 
Wijers, R. A. M. J., 

L51 

Wijers, R. A. M. J., 



1997b, astro-ph/9706060 
E. 1997a, ApJ 485, Lb — 
E. 1997b, ApJ 489, L33 

E., Kulkarni, S. R., & Frail, D. A. 1998, ApJ 497, 288 
D. W. 1986, Quasar Astronomy, Cambridge: Cambridge 



A. 1995, ApJ 453, 583 
Rees, M. J., & Meszaros, 



P. 1997, MNRAS 288, 



Galama, T. J. 1998, istro-ph/9805341 



James E. Rhoads 



15 




2xl0 18 



r (cm) 

Fig. 1. — Dependence of the bulk Lorentz factor T on the burst expansion radius for an isotropic burst and a burst beamed 
into an opening angle ( m — 0.01 radian. Both bursts follow a T oc r~ 3 ' 2 evolution initially, but the beamed burst changes 
its behavior at T m 100 w 1/Cm, beyond which its Lorentz factor decays exponentially with radius. 
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Fig. 2. — Dimensionless peak flux as a function of dimensionless observer-frame time. Points show the results of numeri- 
cal integrations. Da shed l ines show the analytic asymptotic forms from equations |3Cj and |4C|, scaled to fiducial values as 



described in section 3.2.4. The late-time flux density is below the prediction of equation 4(1 by a factor of ~ 0.7. This 



stems from approximations in the exponential regime initial conditions (equations 16 to fL9), which are derived by apply- 
ing the power law regime results beyond their range of strict validity. The solid curve shows an empirical interpolation 
betwee n the early and late-time analytic forms, incorporating a factor 0.7 correction to the late-time asymptotic form (see 
section 3.2.5). 
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Fig. 3. — Dimcnsionlcss frequency at which the spectrum peaks as a function of dimcnsionless observer-frame time. The 
y-axis shows log 10 (^e,mi©)j so that a t^ 2 decay of peak flux density with time appears as a horizontal line. As in figure 0, 
points show the results of numerical integrations; dashed lines show analytic asymptotic forms (here from equations 
and|34|), and the solid line is an empirical interpolation. The late-time asymptotic frequency given by equations j^jis seen 
to be too large by a factor ~ 1.35. This is due to the approximate initial condi tions used to derive equation 34], and a 
correction factor has been applied in deriving the interpolated curve (see section 3.2. 5| ). 
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-5 5 

log lo (t / t eb ) 

Fig. 4. — Sample light curves in dimensionless units. The curves have been derived by combining a broken power law 
spectrum (with electron energy dis tribut ion slope p — 2) with the interpolated, dimensionless forms of the i^.m,©^©) 
and i/® im (t®) relations (see section |3.2.5| for details). The light curves correspond to frequencies 10 8 i/ ffiim ^, 10 1-5 i'® l77l) &, 
and lO^i/ffl.TO.fc. in order of decreasing flux density at the earliest time plotted. All the asymptotic behaviors described 
in equations [5^ are exhibited here, though the additional cooling break (equation |5l]) is omitted. The transition be- 
tween the behaviors for i® <C f©,& and i® ^> i®,6 is rather gradual. The other transition, between i® t® iJ7l (i/®) and 
t® 3> t® >m (i/®), is artificially sharp in these plots because the adopted spectrum has a discontinuous slope at ^©, m - A 
more detailed treatment of the spectrum would smooth out this transition also. 



